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1. Introduction 

Large N gauge theories are very interesting theoretical models. They have simpler properties 
than their finite N counterparts. At the perturbative level only planar diagrams contribute. At 
the non-perturbative level one has factorization, validity of the quenched approximation, stable 
resonances, OZI rule, etc. Furthermore, the connection with string theory is also simpler in the 
large N limit. On the other hand, these theories are expected to be confining and have a rich meson 
spectrum. Thus, it is very atractive to investigate them using lattice gauge theory methods (for a 
recent review see Ref. The problem is that in this context their simple character seems to 
be lost. Indeed, the standard pathway to obtain predictions for these theories is by extrapolating 
the results obtained for finite N. One possible simplification was found by Eguchi and Kawai [||]. 
By examining the Migdal Makeenko equations for Wilson loops, they concluded that, if the 
symmetry of the theory remains unbroken, the large N pure gauge theory on the lattice is volume 
independent. Bringing this idea to the extreme they proposed a reduced single-site formulation 
called the Eguchi-Kawai model. Unfortunately, it was soon realized that in the EK model the 
symmetry is broken at weak coupUng 

Since then, the question has been whether it is possible to make the reduction idea survive 
in the continuum limit. Very soon simple modifications were proposed designed to restore 

the symmetry of the EK model. Recently new proposals have been added to the list 
Alternatively, Narajanan and Neuberger proposed a mixed procedure called partial reduction, 
and suggested that presumably reduction only operates beyond the deconfinement scale. 

Our goal is to test whether the twisted reduction proposal, introduced by the present au- 
thors 1^, is indeed capable of realizing the reduction idea in the continuum limit. This is based on 
the Twisted Eguchi Kawai model (TEK), a single site model whose action is given by 

S = bNY, Ti{z^yU^UyUlul) (1.1) 

where are SU(N) matrices, b is the inverse 't Hooft coupling and z^v are elements of the center 
Z/j?. The choice of z^v is crucial for determining the behaviour of the system in the weak coupling 
limit (large b). Choosing z = 1 gives the original EK model, which has an infinite degeneracy at 
the classical level Although the action is invariant under Z^ symmetry {U^x — z^U^), the 
quantum corrections favour symmetry-breaking minima. 

Other choices of z^v lift the degeneracy of the ground state and the classical vacuum becomes 
invariant under a subgroup of Z^. A particularly elegant choice is termed symmetric twist and given 
by (V > n) 

k 

Zfiv = 4/i = exp{27r/-} (1.2) 

where N = and k is an integer defined modulo L. With this choice, the classical vacuum is 
invariant under a Z^ subgroup of the Z^. In the large L limit (large N) this is enough to secure 
the reduction idea at weak coupling. Perturbative expansion around this vacuum gives rise to 
modified Feynman rules. The propagators become just those of a lattice field theory at finite volume 
L^. This is an important information showing how the N"^ degrees of freedom of the group map 
fully onto effective spatial degrees of freedom. It also suggests what are the dominant finite N 
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corrections to reduction. The remaining Feynman rules introduce momentum dependence at the 
vertices, providing a discretized version of non-commutative field theory. 

Most of the early studies of the TEK model were done for = 1 and showed that the 
symmetry remained unbroken at intermediate couplings. A few years ago it was found [|9|]-[p^]- 



[ ]11[ ] that at intermediate couplings the symmetiy breaks down for N>Nc=lOO. The critical value 
of A'^ depends on k roughly as Nc{k) ~ 90k^. Thus, in Ref. [12] we proposed that the correct large 



A^ Umit has to be taken keeping k/L approximately fixed. Our goal is to test this idea and use the 
reduced model to obtain physical information from the large N continuum gauge theory at infinite 
volume. We are currently involved in this task and the present talk reports part of these results. 

It must be said that we have been able to simulate the model up to A'^=1369 without finding 
any indication of symmetry breaking. Notice that this amounts to effective lattice sizes of 37^ 
which are normally considered large enough in ordinary lattice gauge theory simulations. To try 
to understand the inteiplay between finite A'^ and finite volume in a more detailed form we have 
analysed in certain detail the 3 dimensional version of the theory. Some results from this work, done 



in collaboration with Margarita Garcia Perez, has been reported by her at this conference []13|]. The 
results support our claims that, if the flux k is chosen judiciously, the physical results essentially 
depend on Nl, the product of torus linear size / and group rank. This is indeed consistent with the 
volume independence achieved at infinite N. 

The present work reports another form of evidence in favour of the twisted reduction idea. 
We set ourselves the goal of computing the continuum string tension at large N using the reduced 
model. To make the comparison free of systematics associated with the methodology, we decided to 
use exactly the same procedure to extract the string tension for ordinary SU(N) lattice gauge theory. 



and then extrapolate these results to infinite N. The main results have appeared [ 14 1 recently. Here 
we will show our updated results which include those for SU(4). 

2. Numerical results on the string tension 

The basic data is made up several independent lattice gauge theory simulations for different 
values of N and lattice gauge coupling b and a spatial volume of 32^. Each simulation contains 260 
configurations separated by 100 sweeps formed by 1 heat-bath and 5 overrelaxation updates, after 
discarding 4000 sweeps to reach thermalization. We simulated the SU(N) gauge groups for A'^=3,4, 
5, 6 and 8. For each group we made independent simulations at different values of b (7 for A'^=3 
and 4, and 5 for A^=5,6,8). This huge information can be used to determine the string tension and 
to extrapolate to the continuum limit. 

In addition, we made simulations of the TEK model for A^ = 29^ = S4l, k = 9 and 7 values of 
b. The number of configurations in that case was around 6000 with similar updating procedure and 
separation of sweeps. 

In all cases we computed T x R rectangular Wilson loops W{T,R) computed from smeared 
links, and from them we extracted the Creutz ratios defined as 

y(TR) = -log + '-^'^ + '-^^^(^ - "-^'^ - "-^^ (2 1) 

' ^ ^W{T + 0.5, R-0.5)W{T -0.5, R + 0.5) 

En^ors were computed by jack-knife. Details of the procedure will be given elsewhere [p4|]-[p^]. 
Here we will focus upon the results. 
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Figure 1: compared to the best fit Eq. (2.2) for A?=8 and b=0.3815 



Our goal is to obtain the continuum string tension for all gauge groups. This is a two step 

procedure in which one first obtains the lattice stiing tension as k: = lims j- ycoX{T,R), and then 

one makes a scaling analysis of the results to take the continuum limit. The first part can be done 
by fitting the square Creutz ratio data (R=T) to the formula 



(2.2) 



Different ranges have been examined, and order 1 chi squares per degree of freedom are achieved 
using the data /? = T G [3.5, 8.5]. As an example, we show in Fig. 1 one such fit. 

Once the different values of K{N,b) are extracted from the fits, one must perform a scaling 
analysis. This is done by taking the limit 



. X K(N,b) 
a{N.b)^o a^{N,b) 



(2.3) 



where a{N,b) gives the lattice spacing (in certain units) as a function of the coupling b. Different 
choices of a{N,b) correspond to different renormalization schemes. A good choice makes the 
approach to the a = limit smoother. We have tried several expressions found in the literature for 
a{N,b). In all cases, the approach seems to be linear in with slopes depending on the scheme 
and on the rank of the group A'^. A class of schemes uses perturbation theory and has the advantage 
of expressing the lattice spacing in units of A^. We choose here the choice taken in Ref. [15] 
which has been tested and used for similar purposes as ours. 

Finally, the resulting string tensions o{N) are studied as a function of and extrapolated to 
large N. In Fig. 2 we display the ratio A^/ y^a{N) as a function of 1 /N'^. The data shows a nice 
linear pattern that allows an extrapolation to infinite A'^. The best fit value gives Aj^/ a{°°) = 
0.525(2). The result obtained from the TEK model is A^/ a{oo) = 0.523(5). The agreement 
is very good and provides a nice test that reduction is operative in the continuum within the con- 
finement regime of the theory. The 1/A'^^ dependence displayed by the data matches perfectly with 



the result obtained in Ref. [15|. Our estimate of a{oo), however, disagrees with the one given in 
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Figure 2: A-^/^/a as a funtion of l/N^ 



that reference within statistical errors, although not if the systematic errors given in Ref. [ |15| | are 



taken into account. A new estimate of the large N string tension has appeared recently [|16|]. The 
methodology is also based upon smeared Wilson loops but the volumes and values of used sit in 
between our conventional and fully reduced model. It is interesting to point out that their result is 



in agreement with ours and differs from the best estimate of Ref. []15|]. 

We believe that the main source of systematic error in our estimate of A^/ o{oo) is indeed 
sitting in the scale itself. For example, using other perturbative schemes one gets differences of the 
order of two percent. For that purpose it is better to use a non -perturbative renormalization scheme 
based on the data itself. This can be achieved by inverting the order of the limits R — oo and 
a — > 0. That means that we will take the continuum limit of the Creutz ratios. This is possible 
since Creutz ratios are free from comer and perimeter divergences. Notice, however, that Creutz 
ratios have intrinsic scaling violations due to its mere definition. Thus, in the continuum limit one 
should have 

X{T,R)=a^ {b)F {t,r) + a\b)H{t,r) + ... 



(2.4) 



where t = Ta{b), r = Ra{b) and F{t,r) and H{t,r) are well-defined continuum functions. Indeed, 
F{t,r) is given by 



Fit,r) 



drdt 



a + y{r/t)[^ + ]^]+.. 



(2.5) 



where W{t,r) is the continuum Wilson loop (divergences drop from the formula). The expression 
on the right-hand side is the asymptotic prediction from an effective string theory description of the 
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Wilson loop expectation value. Indeed, the Nambu-Goto theory predicts 

Id/ dlog(ri(iz))\ 

where Tj (iz) is the Dedeking eta function. Indeed, perturbation theory is also compatible with the 
form given in Eq. with a = and 

/ N (jV^ - l)g^ 1 +zatan(z) + atan(l/z)/z 

^''^'^ = iTVz ^ 

The continuum function F{t,r) can be used to fix a physical scale r as follows: 

r2f (r,r) = 1.65 (2.8) 

The choice of 1.65 is conventional and suggested by the analogy between our scheme and the 
so-called Sommer scale. 

Once the scale is fixed, one can use the data to obtain a determination of F{r,t) for each N 
using all the Creutz ratios for all values of 7?, T and b. In particular, our result for r^F{r, r) and A'^=8 

-2 

is given in Fig. 3. The data are very well fitted to by a second degree polynomial in ^ with a fairly 
small quadratic term. Similar behaviour applies for other values of A'^ and for the TEK model. From 
the data we estimate a{oo)r^ = 1.105(10) and 7(1) = 0.272(5). The latter value differs from the 

prediction of Nambu-Goto theory Yng{^) = 0-16 Similarly, we analysed the behaviour of 7(2) 

for z close to one. The data is well described by a pai^ametrization 7(2) = 7(1)(H- 'V 2- )■ ^i^^ 
T = 0.31(6). This value is fairly close to perturbative result 0.39 obtained from 7pr(z)/7f7-(l), and 
differs significantly from the value ~ 2 obtained from Jng- The result might indicate that one has 
to supplement the string contribution with one coming from one gluon exchange. 

3. Conclusions 

The main conclusions of our work are: 

• We have presented a new method to determine the string tension from smeared Creutz ratios. 

• The results scale smoothly to the continuum limit. 

• The continuum string tension extrapolates linearly in 1 /N^ towards the large N Umit. The 
result matches with that obtained from the TEK reduced model. This gives a strong support for the 
validity of continuum reduction. 

• The data satisfy nice scaling properties, that allow the reconstruction of a finite continuum 
function of Wilson loops, called F{r,t). 

• The function gives rise to the definition of a new non-perturbative renormalization scheme, 
and accompanying scale r. 

• The subleading behaviour of F{r,t) differs from the prediction of Nambu-Goto theory. 

A.G-A is supported from Spanish grants FPA2009-08785, FPA2009-09017, CSD2007-00042, 
HEPHACOS S2009/ESP-1473, PITN-GA-2009-238353 (ITN STRONGnet) and CPANCSD2007- 
00042. M.O is supported in part by Grants-in-Aid for Scientific Research from the Ministry of 
Education, Culture, Sports, Science and Technology (No 23540310). 
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(FY2011-12). 
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Figure 3: The continuum function r'Fir^r) for A^=8 
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